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1. Data Analysis Procedures 

A description of the radiative properties of 
two cirrus clouds sampled on 10/28/88 in the 
FIRE cirrus IFO is presented. The clouds are 
characterized in terms of the broadband in- 
frared effective emittance, cloud fractional ab- 
sorptance, shortwave (.3-2.8 ft m) heating rate, 
cloud albedo and vertical velocity. The broad- 
band fluxes used in these calculations were ob- 
tained from measurements made by pyranome- 
ters and pyrgeometers manufactured by Epp- 
ley Laboratories Inc. For a description of these 
radiometers and calibration procedures, see Al- 
brecht and Cox (1976), Albrecht and Cox (1977) 
and Smith Jr. tt ai (1988). In addition, the 
shortwave irradiances were corrected to a hori- 
zontal plane and normalized to the same time by 
taking into account Sabreliner flight information 
(i.e. pitch, roll, heading and angle of attack), as 
well as sun-earth geometry considerations (Ack- 
erman and Cox, 1981; hereafter AC). 

Since only one aircraft was used, broadband 
fluxes at different levels in the cloud were not 
measured simultaneously. As a result, sampling 
errors may occur due to the nonsteady state of 
the cloud field or due to the possibility that the 
flight legs were not flown directly above or be- 
low each other. To minimize these errors and 
to simplify the analysis, the necessary variables 
were averaged and the averages used in the cal- 
culations. The averaging for each flight leg was 
done in the following ways: 

(1) Average the appropriate variables using 
every data point along a flight leg to determine 
a mean value for the cirrus cloud field. 

(2) Average some percentage of the data to 
stratify the cloud into separate regimes (i.e. op- 
tically thicker or thinner than average). 

The downwelling shortwave and longwave ir- 
radiances were used as selection criteria to re- 
move cloud free data encountered along the data 
sampling leg, 

(a) Shortwave Stratification 

Find the lowest 30% of the downwelling short- 


wave irradiance values along each flight leg and 
average all irradiance variables corresponding to 
these points. 

(b) Longwave stratification 

Find the highest 30% of the downwelling long- 
wave irradiance values along each flight leg and 
average all irradiance variables corresponding to 
these points. 

The effective emittance (e*) is the ratio of the 
observed change in irradiance through a cloud 
layer to the change that would have occurred 
had the cloud been black, e*, as derived from 
the downward irradiance may defined as: 

. m _ gg(l) ~ Mk ) 

UJ “ aTg — Ht{1) 

where H( ) refers to the observed downwelling 
infrared irradiance, T is the mean temperature 
along the flight leg, the subscripts T and B re- 
fer to the top and bottom of the cloud layer, 
respectively, and a is the Stefan Boltzman con- 
stant. 

The cloud fractional absorptance (CFA) is de- 
scribed by AC and is defined as: 

CFA T > r Z k zil k Z £» ggj t 

H Jr 

where 

(H I -H t)r = net irradiance at cloud top, 
m -if T)b = net irradiance at cloud base, 
£\- Hsi S net energy gain or loss through cloud 
sides (i.e., dHxfdx + dHy/dy ), 

H ir= incoming irradiance at cloud top. 

In seeking ways to quantitatively express the 
spatial variability of the observed variables, two 
techniques have been applied: a histogram anal- 
ysis and a cumulative variance summary. 

The histograms were created by applying 
thresholds to detrended data and calculating 
the size f in kilometers) of each section that ex- 
ceeded the threshold. These sizes were then 
combined to form the histograms. In Figure 1, 
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an example graph, the +1 SD (standard devi- 
ation) threshold is exceeded between points F 
and G. The distance between F and G is then 
the size used in creating the histogram with the 
threshold of +1 SD. The count of the appropri- 
ate bin size will then be incremented. Similarly, 
the distance between E and H is the size used to 
increment the count of the appropriate bin for 
the +0.0 threshold. For the negative thresholds, 
-0.0 and -1 SD, the values must be less than the 
threshold (instead of being greater than for the 
positive thresholds) for the threshold to be ex- 
ceeded. Thus, the section between points B and 
E, and the section between H and K, exceed the 
-0.0 threshold. Similarly, the section between 
points C and D, and the section between points 
I and J, exceed the -1 SD threshold. Sections 
of the graph at the ends of the data set that 
exceed a threshold (i.e. section to the left of 
point B for the +0.0 threshold) are not used 
because the size is unknown. Each histogram 
contains several different heights to increase the 
total number of sizes. Since the thresholds are 
based on standard deviations and not absolute 
values, the distribution of sizes at each height is 
expected to be similar. 

The cumulative variance graphs were created 
by first using Fast Forrier Transforms to cre- 
ate a variance spectra. The variances were 
then summed in a cumulative fashion from the 
largest scalelengths (wavelengths) to the small- 
est. No detrending was done to the original data 
and no smoothing or averaging was done to the 
spectral points. Also, no points were removed 
due to the effect of the finite length of the file. 
This means that the spectral points at scale- 
lengths larger than about 4000 m (especially 
greater than 8000 m) may only be qualitatively 
correct. The cumulative variance at the smaller 
scalelengths should be correct because a more 
accurate representation of the variance at the 
larger scalelengths should only redistribute the 
energy amongst the larger scalelengths. 

2. Flight Description 

The Sabreliner flight on the morning of 28 Oc- 
tober 1986 was flown over the western shore of 
Lake Michigan. There was a shortwave (pres- 
sure) ridge over the region with a trough ap- 
proaching from the west. 

The cirrus, the Sabreliner flew through, was 
thin and had a banded structure. Cloud base 
of the upper cloud deck varied from about 8.8 
km to 9.2 km with broken undercast beneath. 
Cloud top was at about 11.0 km. 

The Sabreliner flew a racetrack pattern with 
the flight legs being oriented about 30 degrees 


off of the mean wind direction WNW. Because 
of the racetrack pattern, it became apparent 
that the cloud on the NW side had different 
characteristics than the cloud on the SE side 
of the racetrack. Therefore the data from each 
racetrack side were treated separately. “Cloud 
l w , on the SE side of the racetrack pattern, in- 
cluded the flight legs at heights: 11.3 km, 10.0 
km, and 8.8 km. “Cloud 2”, the NW side of 
the racetrack, was sampled at heights: 10.7 km, 
9.4 km, and 8.5 km. Tables la and lb show the 
times, positions and headings of the flight legs 
for each of the designated clouds. 

Table la. Cloud 1 flight leg times, positions 
and headings. 


T*bU 1*. Cloud l flight lag tunas, positions sad heading*. 


( k "») 

Tuns 

(GMT) 

Start 

Lat 

Long 

Tima 

(GMT) 

End 

Lat 

Long 

Heading 

11.3 

15:45:48 

44.29 

-87.89 

15:48:58 

44.45 

-87.41 

120 

10.0 

16:00:56 

44.27 

-87.89 

16:04:31 

44 45 

-87.36 

119 

8.8 

16:16:4? 

44.26 

-87.93 

16:20:41 

44.44 

-87.36 

113 


Table lb. Cloud 2 flight leg times, positions 
and headings. 


Table lb. Cloud 2 Bight leg times, positions sad headings. 


Height 

- (M 

Tim* 

(GMT) 

Start 

Ut 

Long 

Time 

(GMT) 

End 

Let 

Long 

Heading 

10.6 

15:36:22 

44.37 

-87.44 

15:41:54 

44.13 

-88.14 

292 

10.7 

15.53:20 

44.86 

-87.39 

15:57:17 

44.54 

-87 91 

288 

9.4 

1607:45 

44.65 

-87.21 

16:1330 

44.45 

-87.93 

289 

8.5 

16:24:44 

44.52 

•87.36 

16:3034 

44.48 

-88.11 

275 | 


3. Results 

Cloud fractional absorptance values were de- 
duced from pyranometer measurements in the 
.3-2.8 fi m spectral region with a correction 
made for sampling of finite cloud effects in the 
manner described by AC. 

Figure 2 shows a plot of the solar cloud frac- 
tional absorptance vs the infrared emittance for 
the two, two layer cloud samples defined in Ta- 
ble 1. The lower values represent a layer with 
a 1.3 km geometric thickness while the higher 
values represent a 2.5 km thick layer. Values 
derived from both the SW stratification and 
the LW stratification are shown. As we men- 
tioned above, the stratification of the data was 
necessary to eliminate the banded structure of 
cirrus cloud within each data leg. In essence, 
both clouds for both stratifications show very 
similar radiative properties with emittances be- 
tween .48 and .55 for the thicker layers and .28 
and .30 for the thinner layers; similarly the SW 
fractional absorptance varies from .08 to .09 for 
the thicker layers and .05 to .06 for the thinner 
layers. Within the errors of the observations, 
the slopes of the lines connecting these points 
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for the two clouds are essentially the same. 

Layer SW heating rates for these same cloud 
layers are given in tables 2 and 3. The heating 
rates for the two layers are encouragingly similar 
as would be expected from the previous figure; 
for the top 1.3 km layer they vary between .23 
and .29 ° C/hr and for the lower layer between 
.15 and .17 ° C/hr. 

Figure 3 represents the SW reflectance ob- 
served at cloud top for Cloud 2 for three spectral 
bandpasses. Tliis is an example of the structure 
of the observations in the raw data. 

Figures 4 shows the results of the histogram 
analysis of scale size variation on the measured 
vertical velocity. It is notable that most of the 
vertical velocity “cells” have horizontal dimen- 
sions less than 3 km. Figure 5 shows that the 
downward visible irradiance shows essentially 
the same structure as the vertical velocity. 

Figure 6 shows the cumulative variance anal- 
ysis for the albedo at the top of cloud 2. The 
continuous line is a fit to the plotted data points. 
This graphs may be interpreted to show the per 
cent variance explained by scale lengths larger 
than a given value. For example, this figure sug- 
gests that 85% of the variance arises from fea- 
tures whose horizontal scale is greater than 10 
km. Figure 7 illustrates the cumulative variance 
analysis for vertical velocity applied to each of 
the flight levels in the sample of cloud 2. The 
solid curve represents the 10.6 km height, the 
dashed curve - 9.4 km, and the dotted curve - 
8.5 km. Notable is the narrowing of the scale 
length with increasing height. Fifty percent of 
the variance is explained by features larger than 
13 km at the 10.6 km height, larger than 22 km 
at the 9.4 km level and larger than 50 km at the 
8.5 km level. 
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TkbU 2. Layer hniik| rttw (.3*2.1 pm) in degreee/hr. for Cloud I 


Z * 1 1.3 Km 

MEAN 

1R STRATIFICATION 

SW STRATIFICATION 


• 

.233 

.261 

Z m 10.0 Km 

• 

(ji«) 

(.309) 


•ITI 

■ IM 

.173 

tm 8.8 Km 

(■»>) 

(M») 

(MI) 


Valuta in 

parentheses uncorrected for cloud inhomogeneitie* 1 


'indicate* heating rate indeterminate 


Table 3. Layer heating ratce (.3-2.8 pm) in degrees/hr. for Cloud 2 


Z a* 10.7 Km 

MEAN 

IR STRATIFICATION 

SW STRATIFICATION 


• 

.289 

.266 

Zm 9 A Km 

• 

• 

(.431) 


.171 

.166 

.148 

Z « 8.5 Km 

(•338) 

(»*) 

(.158) 


Value* in 

p area theses uncorrected for cloud inhomogcnciitce j 



Figure I. Example plot illustrating procedure 
for generation of size histogram plots (figs. 
4 and 5). 
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FIRE CIRRUS IFO 26 OCT 86 FLT. 12 



Figure 2. .3-2.8 p m CFA va. effective emittance 
for Cloud 1 and Cloud 2. 


FIRE CIRRUS IFO 28 OCT 86 FLT. I 2 
ALBEDO (10.7 km) 



Figure 3. Albedo for Cirrus Cloud 2. 


FIRE Cirrus IFO 28 OCT 86 Fit, 12 
Vertical Velocity 



Figure 4. Histogram illustrating observed verti- 
cal velocity scale sizes (SD is standard devi- 
ation from detrended mean). 



Size (km) 

Figure 5. Same as Figure 4 but for downwelling 
shortwave irradiance. 


FIRE Cirrus IFO 25 OCT 86 Fit. 1 2 



Scoiefength (m) 


Figure 6. Cumulative variance of .3-2.8 ^ m 
albedo for Cloud 2. 


FIRE Cirrus IFO 28 OCT 86 FH. 1 2 



Scoietength (m) 


Figure 7. Cumulative variance of vertical veloc- 
ity at different heights in Cloud 2. 









